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ABSTRACT 

Gimballess inertial reference equipment might provide six possible 
advantages over conventional stable member systems. These are savings 
in power, size, weight, and expense, flexibility of packaging, and ease 
of maintenance. The theoretical problems involved in gimballess equip- 
ment are the question of what inertial instruments measure, and how 
best to use the information they generate. The engineering problems 
concern the difficulty of torquing gyroscopes with accuracies approaching 
one part in a hundred million, and the design of inertial instruments to 
withstand the high angular rates they will experience in this application. 
The inertial package to be considered contains three delta-modulated inte- 
grating gyros and three delta-modulated integrating accelerometers. 

Delta modulation is a pulse modulation system which delivers 
output pulses, representing increments of the input or one of its 
integrals or derivatives, in synchronism with a clock. The analysis 
of delta -modulated loops with either a two- or three -level relay is 
a straightforward matter using z-transforms and modified step-by-step 
transient methods. Application of delta modulation to integrating gyros 
and accelerometers leads to instruments whose output pulses represent 
increments of rotation and velocity. 

Transformation of a vector from one Cartesian coordinate frame 
to another may be represented in three essentially different ways: 
through matrix-tensor operations, four-parameter techniques, and 
vector algebra. Comparison of these to determine which is most 
suitable for machine computation, both from the standpoint of generating 
the transformation parameters and of accomplishing the transformation, 
shows that a matrix transformation, with straightforward generation of 
direction cosines, is best. The most attractive choice of machine for 
the computation is the digital differential analyzer, since the angular and 
linear motion information is in incremental form. High-order integration 
rules with larger angular increments cannot be used to reduce computer 
speed while maintaining accuracy because of an irreducible uncertainty 
in the digital angular information, 



The transformation computer may be analyzed using frequency domain 
techniques. With these techniques, a signal flow graph representing the 
computer can be developed, which gives great insight into the dominant 
mode of operation of the computer, and which provides an estimate of 
truncation error. Quantization may be treated statistically. 

This investigation has considered one way of instrumenting gimballess 
inertial reference equipment in detail, and has provided a survey of some 
of the possibilities in this field. It has shown the theoretical feasibility 
of such equipment, and has also indicated that the major development 
effort required to permit practical realization of such equipment with 
a specified accuracy must be concerned with inertial instruments designed 
especially for this application. 
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CHAPTER 1 



INTRODUCTION 



l.A Background 

The problem of arriving at a destination has been important to man 
from the beginning. At first, getting there was merely a matter of 
following a line of sight to the destination. Later, landmarks whose 
position with respect to the destination were known were used as 
guideposts in making journeys out of sight of the destination. The stars, 
which can be used out of sight of land, represent an important class of 
these reference points. Celestial bodies can be observed only on clear 
nights, and of course, other landmarks, by and large, must be visible 
to be useful. With the invention of radio, and later radar, man-made 
radiation could be used to distinguish landmarks during periods of 
reduced visibility. In modern military and transport applications many 
situations exist in which terrestrial and celestial reference; points 
are either unavailable or unreliable as a basis for guidance. In such 
cases there is a need for navigation systems which do not depend on 
any external information. Inertial navigation systems, which deduce 
the position of a vehicle from on-board measurements, are systems 
of this type. The measured acceleration is integrated twice with respect 
to time to give displacement for an initial position. 

The majority of inertial navigation systems today use two or three 
single-axis accelerometers mounted on a stable platform to measure 
acceleration. The platform is either kept non -rotating with respect 
to inertial space or precessed so that the input axes of the accelerometers 
are kept coincident with some slowly rotating set of coordinate axes, 
such as local geographic coordinates. The stabilized platform is 



articulated with respect to its support by a system of gimbals. The 
gyros mounted on the platform detect the angular motion of the platform, 
and the gimbals are torqued in response to gyro signals to drive the 
platform back to its null position. 

Another configuration for inertial navigation systems has the inertial 
instruments fixed to the vehicle and hence is gimballess. Single axis 
accelerometers measure vehicle acceleration and resolve it into vehicle 
coordinates. The guidance equations may be solved in body coordinates, 
in which case gyro information about vehicle rotations is used directly 
in the guidance equations; or they may be solved in an inertial coordinate 
frame, in which case the gyro information is used to generate a coordinate 
transformation which operates on the accelerometer data to give acceleration 
resolved into inertial axes. 

There are six possible advantages available in gimballess guidance 
systems. With the removal of gimbals, weight, space, and power 
required are reduced. Further, the not inconsiderable expense of 
gimbals is eliminated. Since the inertial package in a gimballed system 
must assume practically a spherical shape, it is necessary to find a 
spherical space in the vehicle in which to put the guidance system. 

Without gimbals to dictate packaging, much more flexibility in packaging 

is available. Finally, assembly and maintainence of the inertial package 

is easier without the encumbrance of gimbals. 

l.B The Problem of Gimballess Inertial Reference Equipment 

1. B. 1 Theoret ical As pects 

The theoretical problems of gimballess inertial navigation are 
two: the measurement of acceleration in body fixed coordinates and 
the determination of the orientation of the body fixed coordiante frame 
with respect to the inertially fixed coordinate frame from gyro information. 

The three single axis accelerometers measure inertial acceleration 
and resolve it into components along body-fixed axes. Coriolis effects 
are not inherently present in the measured acceleration components; 
they arise only in the intepretation of these components in non-inertial 
reference frames. If such frames are used in the guidance equations, 
the Coriolis terms must be computed separately using angular information 
from the gyros. If the guidance equations are set in inertial coordinates. 
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the three measured components need only be transformed into the inertial 
coordinate system. 

The generation of a vector transformation from gyro information is 
a straightforward process. The differential equation governing the 
transformation parameters may be derived and instrumented in any 
convenient way. 

1. B. 2 Engineering Aspects 

The engineering problems involved in gimballess applications of 
inertial guidance are many, but the two most important involve the 
environment to which the inertial instruments will be subjected and the 
speed of computation required to keep the computed orientation of the 
body -fixed coordinate frame close to its actual orientation. 

While inertial instruments mounted on a stable platform are isolated 
from the angular motions of the vehicle carrying them, they do experience 
all the linear motion and most of the vibrations of the vehicle. Body- 
mounted instruments, in contrast, experience both the angular and linear 
motions, but they may be effectively shock mounted to minimize the 
vibrations present. Although much effort has been expended in developing 
instruments whose performance does not deteriorate under acceleration, 
practically none has been devoted to designs which might experience large 
angular velocities. For this reason no instruments are available at 
present which can be used in a gimballess application. 

In gimballess inertial reference equipment, the gyros are used to 
measure angular orientation of the body-fixed coordinate frame with 
respect to some reference. This may be done either by using a pick-off 
arrangement with free gyros, or by torquing single-degree-of-freedom 
gyros to make them follow body rotations. If free gyros are used, the 
pick -offs must be accurate to about five thousandths of a degree. The 
required accuracy of torquing for single-degree-of-freedom gyros may 

be shown by the following considerations. If the gyro is to have a drift 

_ 8 

rate of one meru (about 7 .10 radians per second), the uncertainty in 
applied torque must not exceed this level. The maximum angular velocity 
which a typical vehicle might experience is about one radian per second, 
so that the maximum torque available must be able to precess the gyro 
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at that rate. Comparison of these two numbers shows that for gimballess 
applications, and the above specifications, the torquing must be accurate 
to about seven parts in a hundred million. For higher acceptable drift 
rates or lower maximum velocities,, this accuracy requirement is less 
stringent. 

The orientation of the body-fixed coordinate frame must be known 
to a high degree of accuracy (about five thousandths of a degree); thus 
at reasonable angular velocities very high computer speeds are required 
to keep up with the vehicle rotation. If a single adder is to be used, 
add times on the order of five microseconds per word are required. 

This is near the limit of available performance. 

1. C Scope of the Investigation 

A gimballed inertial reference system delivers acceleration information 
in stabilized inertial coordinates. This investigation will consider an 
inertial reference system with single axis integrating accelerometers 
and single-degree-of-freedom integrating gyros fixed to a vehicle. All 
the instruments will be pulse torqued. The gyro information is to be 
used to generate a coordinate transformation which will operate on 
accelerometer information to give inertially stabilized data. Thus it 
is proposed that the gimballed inertial package be replaced by a body- 
mounted inertial package and a transformation computer. 

This investigation assumes that the guidance equations will be 
solved in inertial coordinates; no examination is made of the possibility 
of doing all or part of the computation in body-fixed coordinates. Further, 
no consideration is taken of the feedback through the guidance computer 
and vehicle autopilot. 

There are no studies available which compare the attractiveness 
of various coordinate frames available for computation in gimballess 
guidance. For this reason, the present investigation of the performance 
and design of a system which provides inertially stabilized information 
from a gimballess inertial package is a useful one. Further, since the 
action of the inertial package and transformation computer is much 
faster than the feedback present, the equipment may profitably be examined 
as an open loop device. 
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Since the inertial package consists of three delta-modulated single- 
axis integrating accelerometers and three delta-modulated single-degree - 
of-freedom integrating gyros, the theory of delta modulation and its 
application to inertial instruments will be considered in some detail. 
Methods for deteiminLng the zero-input Dehavior of deita-moduiated loops 
will be reviewed, and the response of such ..oops to impulses and steps 
will be found. Criteria for the choice of loop parameters will be set 
forth. 

The available methods of transforming a vector from one Cartesian 
coordinate frame to another will be examined. There are essentially 
three ways of viewing the transformation: the matrix of direction cosines, 
the four parameter representations, and the vector representation. The 
properties of each will be considered and the differential equations 
governing the transformation parameters will be derived. Finally, the 
transformations will be compared to determine their suitability for 
machine computation, considering the form of instrument output. 

The transformation computer will be designed to instrument the 
transformation found to be best for machine computation. The possibility 
of accumulating a number of increments from the gyros and instrumenting 
more complicated algorithms to permit the use of a computer with a 
slower clock rate and reduction of error will be examined. With the 
design of the computer fixed, a detailed analysis will be carried out 
to determine the dominant mode of the computer and estimates of 
truncation and quantization errors. 

Finally a few typical values will be used to illustrate the method 
of determining the basic system parameters and at the same time to 
provide a feel for the order of magnitude of the numbers involved. 

1. D Previous Resuits 

Delta- modulation was originated as a transmission system (Refs. 

2, 5, 6, 14) with little attention being paid to the feedback aspect of 
the transmission loop. Later, with interest developing in contro. 
systems using relays, control theory was applied to deita-moduiated 
loops (Refs. 1, 3, 4, 6-13, 15-19) under the names bang-bang, on-off, 
and contactor control systems. The zero-input behavior of these loops 
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received most of the attention, with phase plane methods playing an 
important part (Refs. 8-11). In 1950 a method for determining the exact 
zero input behavior of a loop of any order by means of difference equations 
was presented by Torng and Meserve (Ref. 19) and the z-transform 
equivalent of this method was suggested by Bergen in his discussion of 
their paper (Ref. 1). Schwartz elaborated Bergen's method and applied 
it to the determination of the loop's response to step inputs (Ref. 15). 
Beyond this, no work has been done previously in determining the response 
of the loop to step inputs . 

The application of delta modulation to inertial instruments first 
took place in 1957, and development work has continued so that these 
instruments are now a very attractive choice for any inertial system 
which requires instrument outputs in digital form (Refs. 4, 7, 8, 17, 18). 

Coordinate transformations have received only cursory examination 
in connection with gimballess inertial reference equipment. Although 
some reports consider the possibility of using transformations other 
than direction cosines, they do not examine them in any detail. Of 
course, the literature on coordinate transformations is available in 
almost any text or reference on classical mechanics or mathematical 
physics (Refs. 30-39). 

The analysis of numerical processes has a fairly long history, 
and most techniques are easily applied to digital computers (Refs. 

41, 44, 46, 48, 51-53). However the application of frequency domain 
methods to the analysis of computers is quite rare (Refs. 49-51, 55, 61), 
and thus most analysis to date is quite cumbersome. No attempt has 
yet been made to apply frequency domain analysis to digital computers 
in which the vairable of integration has not been time. 

A number of private companies have worked in the field of 
gimbaLless guidance, both on their own initiative and under contract. 

Their studies have mostly provided a brief survey of the possibilities 
in this field, but no detailed analyses. Very little of this work is 
available in the open literature (Ref. 66). 
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CHAPTER II 



DELTA -MODULATED LOOPS 



2. A Introduction 

Delta modulation is a puise modulation system patented in France 
in 1946 (Ref. 6) and first described in the literature by de Jager (Ref. 5) 
and others (Ref. 14) in 1952 as a system for speech transmission. In 
1957, a similar system was developed giving digital outputs directly 
from instruments and interest in such instruments has remained active, 
(ftefs. 4, 7, 8, 12, 17, 18) 

2.B Description of the Loop 

The input is translated into a sequence of pulses, each having the 
same value, but with either polarity being possible. The pulses are 
obtained by means of the negative-feedback loop shown in Fig. 2-1. One 
or the other of the iinear plant functions may be unity. When the loop is 
used for transmission, the iinear plant in the forward path is usually 
unity, while in instrumentation applications the feedback plant usually 
contains a zero-order hold. In any case, an integration occurs somewhere 
in the linear plant to assure zero average error (Ref. 3, p. 53). The 
relay may have either two levels or three levels. The mathematical 
model for the sampler is an impulse modulator synchronized with a clock. 

Each impulse appearing at the output represents a unit change in the 
input or one of its integrals or derivatives. This representation of a 
signal in terms of increments leads to the name "deLta modulation". If 
the integration appears in the forward path, each output pulse represents 
a change in the integral of the input. If the integration appears in the 
feedback path, each pulse represents a change in the input itself. Two 
integrations in the feedback path cause each pulse to represent a change 
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Fig. 2-1 Delta- Modulated Loops 
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in the slope of the input. 

2.C Analysis for Zero Input 

With no input, a loop with a two-level relay will experience a self- 
excited oscillation or limit cycle, and if the zero level is narrow enough, 
a loop with a three-level relay may also have a limit cycle. From physical 
considerations, it is clear that the periods of these limit cycles are 
integral multiples of the sampling. period. Further, if the linear plant 
contains an integration, the number of positive pulses in a cycle will 
equal the number of negative pulses (Ref. 3, p. 53). 

2.C.1 Moding 

Each limit cycle exhibited by delta-modulated loops produces a 
characteristic sequence of output pulses. These sequences are referred 
to as the modes of the system and are very important in the analysis of 
input-output relationships. The modes of a system are designated by 
the output sequence they produce. For example, if the limit cycle in a 
loop with a two-level relay produces a single positive pulse followed 
by a single negative pulse, the mode designation is 1-1. Similarly, a 
mode producing two positive pulses followed by two negative pulses is 
designated 2-2. If the loop contains a three -level relay and has a limit 
cycle producing two positive pulses followed by two negative pulses and 
one zero pulse, the mode is termed 2 -0-2-1. The zero for the second 
digit indicates that no zero pulse occurs between the positive and negative 
pulses. Note that the mode designations always begin with the positive 
pulses. 

The concept of moding is applicable only to systems with no input, 
since the limit cycle is disturbed by inputs. The description of outputs 
in response to inputs is discussed in Section 2.D. 

2. C. 2 Describing Function Analysis 

In his famous paper(Ref. 58) Kochenburger presented a method of 
steady-state sinusoidal analysis applicable to non-linear systems in 
which non-linear elements are represented as linear gains. These 
gains are determined as functions of the amplitude and frequency of 
the inputs to the non-linear element. Chow (Ref. 3) and Russell (Ref. 13) 
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first used the describing function to determine limit cycles in delta- 
modulated loops. Their interest was in adjusting the width of the zero 
level in a three-level relay to eliminate limit cycles, but their analysis 
is applicable to any sampled-data loop containing any non-linearity. 

Their method of analysis is the same as the standard describing 
function analysis applied to continuous time systems with the addition 
of a variable delay due to sampling. This variable delay, 6, arises because 
the output pulses can occur only at sampling instants, even though the 
relay had switched earlier. This delay, which ranges from zero to a 
full sampling period, T , is easily understood from an examination of 
Fig. 2-2. The marks on the time axis indicate clock pulses. In each 




Fig. 2-2 Time Delay due to Sampling 

case, the negative output pulse occurs with the second clock pqlse shown 

\ 

and 6 ranges from almost an entire sampling period in Fig. 2 -2(a) to 

practically zero in Fig. 2 -2(c). The maximum phase shift due to sampling, 

<j) , is the product of frequency and sampling period, and is given by 
s 

T 

0 S = W T S = 2? r T = M'Z (2-1) 

s 

where T is the period corresponding to the input angular frequency o>. 

Standard describing function analysis (Ref. 63, Ch. 10, pp. 559-580) 

is carried out by graphical solution of the equation 

1 + G(j«)N(A, w) = 0 (2-2) 

G(jw),the transfer function for the linear part of the system, is plotted 

1 

on a Nyquist diagram with - the negative of the reciprocal of 

the describing function for the non-linear element. Intersection of the 
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two curves indicates a root of Eq. (2-2) and predicts a limit cycle. In 
the analysis of delta-modulated loops the phase shift due to sampling 
must be added to the plot; this is most conveniently done by adding 0 

s 

as a phase lag at frequencies on the G(jw) locus for which a limit cycle 
is possible, that is, at submultipies of the sampling frequency. In this 
application, a log polar plot is most convenient for making the Nyquist 
plot. It is termed a "moding plot". ^ 

In Fig. 2-3, a moding plot for G(jw) = j w (j WT + is shown. The 

negative reciprocal describing function for a two-level relay covering 

the whole negative real axis is shown. As can be seen from this plot, 

limit cycles are predicted for angular frequencies of o> g /2, w g / 4, and 

to /6. At to /8 the sampling phase is at mt>st tt /4, which is not sufficient 
s s 

to reach the negative real axis. It is clear that limit cycles will not 
exist at lower frequencies; thus the modes predicted are 1-1, 2-2, and 
3-3. For a three-level relay, the describing function reaches a maximum 
when the amplitude of the relay input is n/ 2 1 times the zero level half-width, 
A c » Thus the reciprocal of the describing function has a minimum, and 
the zero level half-width may be adjusted so that no limit cycle is 
predicted. An example of this calculation is shown in Appendix B. 

2. C. 3 Detailed Analysis 

The describing function is only an approximate tool, albeit a powerful 
one. For this reason attempts were made to achieve exact answers. 

The step-by-step transient method was discarded as being too cumbersome. 
Phase-plane methods received more attention, and although limited to 
second order systems became quite elegant (Refs. 9-11). Torng and 
Meserve (Ref. 19) presented a method for determining the exact output 
of the linear plant of a delta-modulated loop at sampling instants using 
discrete orthogonal functions. Bergen (Ref. 1), in his discussion of 
their paper, outlined a similar method using z-transform techniques. 

His method allows the possibility of determining no-transient initial 
conditions for the limit cycle. By manipulation, the state of the linear 
plant at each sampling instant may also be determined. 

Each method assumes a mode at the relay and applies it as an input 
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to the linear plant. The output of the linear piant resulting from this 
input is examined at the sampling instants to see if the signal has the 
proper polarity to sustain the postulated mode. Thus if the 1-1 mode 
be possible, the linear piant output at sampling instants must be a 
repeating sequence of two numbers: the first positive, the second negative. 
Torng and Meserve do this using discrete orthogonal functions in the time 
domain. Bergen operates in the frequency domain with z -transforms, and 
introduces initial conditions in this way. 

The modes which must be investigated are limited, and may be 
determined by a describing function analysis as carried out above ( Sec » 

2. C. 2). Torng and Meserve list the modes which must be investigated 
for limit cycles with periods of two, four, six, and eight sampling periods 
(Ref. 19, Table II, p. 303). 

For a loop with a three -level relay, choice of the zero-level width 

so that a describing function analysis predicts no limit cycle may be 

made. However, for fast systems (T / T > 1), this choice of zero-level 

s 

width may not eliminate low energy limit cycles. To eliminate these, 
a transient analysis is the easiest. The linear plant output is assumed 
to be just at the edge of the zero level of the relay. Thus one feedback 
pulse drives the linear plant back toward null. For no limit cycle to 
occur, the limit of the linear plant motion must not reach the outer 
edge of the zero level, even for very long times. The result of this 
analysis for any linear plant containing an integration is that the zero 
level width should be chosen greater than the product of the sampling 
period and the steady state velocity of the plant under the action of the 
relay alone. 

2.C.4 Compensation 

As Schwartz (Ref. 15) and Sliwkowski (Ref. 16) point out, in limit 
cycling loops inputs may not be sensed until one full period of the 
prevailing mode has elapsed. Another view of this problem is that in 
limit cycling, the loop is storing information; the amplitude of the limit 
cycle, indicated by the moding, is a measure of the information storage. 
This information storage should be kept to a minimum if the system 
is to be used to best advantage. In a loop with a three-level relay the 
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width of the zero level is a measure of stored information, and for slow 

systems (T /t < .1) may become quite large. For this reason it may be 
s 

desirable to compensate the system so that all but the 1-1 mode are 

eliminated, or so that the zero level of a three -level relay is reduced 

to a reasonable size. There are three methods of compensation which 

suggest themselves: dc analog, ac analog, and discrete or digital. 

Compensation is best viewed in terms of the moding plot (Fig. 2-3). 

To eliminate the 2-2 mode, more than 90° phase margin must be provided 

for G(jw) at w /4. Similarly to eliminate the 3-3 mode 60° of phase 
s 

margin must be provided at w /6 and so on. With these conditions met, 

s 

the extended G(jw) locus will not intersect the describing function except 
at the 1-1 mode. For a three -level relay, the zero level width may be 
adjusted so that the reciprocal describing function reaches down to just 
above the 1-1 level. 

Dc analog compensation methods are covered in detail in most 
elementary automatic control textbooks (Ref. 57, Vol. I, Ch. 12; Vol II, 

Ch. 4; Ref. 64, Ch. 6; for example). Adjusting the moding plot to achieve 
the objectives listed above is a straightforward task. Ac analog compensation 
is also treated in many texts (Ref. 63 Ch. 6; Ref. 64 Ch. 6; for example). 

The desirability of this form of compensation stems from the fact that 
the relay is often a phase-sensitive switch operating on information in 
a supressed-carrier wave; thus, if dc compensation were used, the 
signal would have to be demodulated. However, the signal frequencies 
used in practice are very close to the carrier frequencies, and de- 
modulation is very difficult under these conditions. For this reason, ac 
compensation is quite attractive. It should be noted that the use of low- 
pass to band-pass transformations do not provide sufficiently accurate 
compensation foi’ this application. 

Discrete or digital compensation has received little attention in the 
literature. The advantage of this form of compensation is that there is 
no problem of drift as compared to analog compensation, and non-linear 
or logical compensation is possible. Most books on sampled -data systems 
include sections on the techniques of linear discrete compensation and 
Tou (Ref. 62, pp. 444 ff) lists three methods for realizing linear discrete 
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transfer functions: computers, delay lines, and pulsed RC circuits. 

The point in the delta -modulated loop at which the signal is discrete 
is after the sampler; it is here that the discrete compensator must be 
inserted. If a pulsed RC circuit is used, a second sampler must follow 
it. 

Linear discrete compensation is best viewed in the G(z) plane. The 
locus of the linear plant has extensions similar to those in the moding 
plot, because of sampling delays. The compensation must be chosen so 
that the extended G (s) locus does not intersect the negative reciprocal 
describing function locus. 

Logical discrete compensation, using branching operations which 
influence the relay on the basis of past output pulses, has been tried only 
once (Ref. 20) with little success. It would seem that there are great 
possibilities for this type of compensation. 

2.D Input -Output Relationships 

As explained above, each output pulse represents a unit change in 
the input or one of its integrals or derivatives. Thus with a finite input, 
the algebraic sum of the output pulses should average to the input or one 
of its integrals or derivatives or, in the case of zero input, to zero. 

The system overloads when the rate of change of the input (or one of 
its integrals or derivatives) exceeds that value which requires all the 
output pulses to be of one sign (Ref. 2, p. 145). Therefore a convenient* 
way of expressing the input is as a fraction of the maximum permissible 
signal. 

Since the position and number of integrators in the loop change input- 
output relationships, a general analysis is very cumbersome. For this 
reason, the analysis discussed here considers a loop with one integrator 
in it. The linear plant in the feedback path is moved to the forward path, 
with the exception of zero -order hold, if one is present, which is left 
in the feedback path. This requires the reciprocal of the moved plant 
to be placed in front of the loop. Fig. 2-4 shows the rearranged loop. 

The dotted lines around the hold indicate that it may or may not be 
present. G!>(s) is the part of the linear plant moved out of the feedback 
path. The analysis which follows ignores the l/GJ, (s) in front of the 
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Fig. 2-4 Delta-Modulated Loop, Rearranged 
2.D. 1 Transient Response 

The following analyses are exact determinations of the inpulse and 
step responses of delta modulated loops. With the exception of Schwartz's 
work (Ref. 15) on step inputs for loops with two -level relays, they have 
not been previously described. 

2.D.l.a Impulse Response 

The height of an impulse (infinity) exceeds the capacity of the system; 
the best output makes up the integral of the impulse a step at a time to 
within the quantization error of the loop. Clearly if the value of the 
impulse is too low, no output will appear, except for the limit cycle. 

The analysis of the loop's impulse response is best done by a step- 
by-step transient method. The applied impulse drives the linear plant 
Some distance from null; calculation of this excursion in the presence 
of the feedback signal is a simple task. The linear plant output at 
sampling instants is examined for the time at which the feedback signal 
changes sign. If the loop has a two-level relay, the state of the linear 
plant is compared with the set of initial conditions corresponding to 
each mode to determine into which mode the system will settle. Usually 
this will be the highest mode. The output of the linear plant plant in a 
loop containing a three-level relay is examined to see if it will cause 
the relay to have an output of opposite sign from the input. If not, the 
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computation is complete; if so, it must be carried further. 

2.D.l.b Step Response 

A step input of magnitude D, where D is the drive level of the relay, 
will saturate the system. For this reason the magnitudes of step inputs 
are given as a fraction, r, of D. 

For a periodic sequence to appear at the output, the input must be 
expressed as a rational fraction of the saturating input. This is not a 
serious limitation, since any number less than one may be well approximated 
by a rational number with a small denominator 

Suppose the input is a rational fraction, say U/U, of the saturating 
signal. Then the output will be a sequence of NW pulses. For a loop 
with a two-level relay, there will be ]Nf + = (W+U)/2 positive pulses and 
N = N(W-U)/2 negative pulses. This may be shown by a short computation. 
The relationship between U, W, N + , and N is 

N + -N~ _ U 

N + +N" W (2-3) 

Rearrangement of this ratio yields 

N + _ W+U 

N W-U (2-4) 

A loop with a three-level relay may have N° zero pulses. In this case, 
the required relation between U, W, N + , N , and N° is 

N + -N~ = JJ 

N + +N~+N° W (2-5) 

If U is positive, and there are no negative pulses, then 

N + U 

N° W-U (2-6) 

With ratios, of positive, negative, and zero pulses determined, it only 
remains to determine the organization of these pulses within the sequence. 
This may be done by determining bounds on the number of consecutive 
pulses of one sign under various conditions. Different bounds are of 
interest depending on whether the loop has a two-level or a three-level 
relay. 
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For a loop with a two -level relay, the bounds that are important are 
the maximum number of pulses of one sign that are possible and the 
minimum number of pulses of one sign that are possible following a 
given number of pulses of the opposite sign. 

The upper bound is found by observing the greatest possible excursion 
from null and counting the number' Of clock periods required to drive the 
linear plant back across null. The maximum excursion from null occurs 
when the linear plant output travels away from null at maximum rate for 
a full sampling period. 

The lower bound is found by assuming minimum initial conditions. 

The minimum displacement is clearly zero plus, and the minimum 
derivatives result from maximum velocity at the previous switching 
time a given number of sampling periods earlier. Under these initial 
conditions, the number of sampling periods required for the system to 
cross null is found. 

The bounds to be determined in the case of a loop with a three-level 
relay are the maximum number of pulses of the same sign as the input 
and the maximum number of opposite sign. 

Determination of the maximum number of pulses of the same sign 
as the input proceeds from the assumption that the linear plant has 
travelled at maximum rate under the influence of the input alone for a 
full sampling period. At this point, positive pulses begin to drive the 
plant back to null. The pulses required to bring the linear plant output 
within the zero level are counted. 

The possibility of pulses of a polarity opposite to the input may be 
examined by a similar analysis. It is assumed that the linear plant is 
driven across the zero level under the influence of the feedback signal 
and input for a full sampling period starting from just outside the zero 
level. The feedback signal switches off, and the system output is 
examined at successive sampling times to see whether it gets outside 
the zero level on the other side. 

Examples of the above computations are shown in Appendices A and B. 

With these bounds set, the sequences having the required average 
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may be written down. These sequences are designated in a way similar 
to that used for modes. For example, a sequence from a loop with a 
two -level relay consisting of three positive pulses followed by a negative 
pulse is designated 3, 1; one having four positive pulses, one negative 
pulse, two positive pulses, and another negative pulse is designated 
4, 1, 2, 1. A loop with a three -level relay having an output sequence of 
a positive pulse followed by a zero pulse would have that sequence de- 
signated 1, 1, 0, 0. Nothe that each of these sequences has an average 
of one -half. The permissible sequences, which are few because of the 
elimination process just described, may be tested by Bergen's method 
to see if they can actually occur. The sequence, with the input, is assumed 
to drive the linear plant, and the signs of the linear plant output are 
examined at sampling times to see if the postulated sequence is sustained. 

With inputs smaller than about one -tenth of the saturating signal, 
the period of the output sequence becomes rather too long for the analysis 
above to handle conveniently, and a modified step-by-step transient 
procedure is preferable. 

In a loop with a two -level relay and a small input, the prevailing mode 
will continue for a number of cycles. The feedback signal has zero average 
over a full period of the mode, and its effect may be ignored. This is 
perhaps best visualized in phase space with the mode trajectory moving 
slowly through the space under the influence of the input. When the 
position coordinate of one of the sampling points of the trajectory moves 
to a new region of phase space (from positive to negative, for example), 
the prevailing mode will change. The mode changes either by dropping 
a pulse or by adding a pulse. In any case, the total state of the linear 
plant after the mode changes may be found and compared with the initial 
states of the other modes. Thus the new mode of the system is found. 

The trajectory of the new mode is assumed to move under the influence of 
the input, with a new initial position and an initial velocity which is the 
velocity of the plant relative to the new prevailing mode. This process 
is continued until the entire sequence is determined. 

For the loop with a three -level relay, the appropriate calculation 
is a step-by-step transient response. Assume an initial condition, and 
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follow the plant through until the sequence is determined. In almost 
all cases only one positive output pulse will occur between zero pulses. 
2.D.1.C Effects of Imperfect Integration 

Thus far, it has been assumed that the integrator has been perfect. 

In some cases, the integrator is in reality a first order lag with a. long 
time constant. This may be viewed as a small elastic restraint. With 
an elastic restraint present, there is a limit below which inputs cannot 
be detected. This is also best seen in terms of the mode trajectory in 
phase space. If the steady-state response of the linear plant to a small 
input is insufficient to move the trajectory far enough to cause the mode 
to change, no input will be detected. Of course, with a three-level relay, 
this critical value of linear plant output is the edge of the zero level. 
2.D.l.d Ramps and Higher Order Inputs 

Ramps and higher order inputs (e.g. parabolic) cannot be followed 
by the loop over infinite time sipce the amplitude of such a signal increases 
without limit. A sawtooth wave may be analyzed in terms of the frequency 
response described below. For ramps which rise to a certain value and 
then remain constant, a step-by-step transient computation seems to be 
the most convenient. The bounds on number of pulses are useful in this 
computation. 

2.D. 2 Frequency Response 

Since the frequency response method of servomechanism design is 
familiar to most engineers, the following modification of the describing 
function is presented. 

A sinusoidal signal is assumed at the relay. The standard describing 
function is used to represent the relay, assuming enough low pass filtering 
in the loop to make the describing function valid. The fundamental of the 
relay is passed through the sampler and hold, giving the feedback signal. 
With the input to the relay, which is the output of the linear plant, known, 
the error signal can be determined. The input to the loop is the sum of 
the error and feedback signals. With the relay output taken as the loop 
output, the describing function for the entire loop is determined as the 
ratio of relay output to loop input. 

Three comments must be made about the above analysis. First, 
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the output of the loop has been taken as the relay output, with no limit 
cycle present. Thus, this analysis, for frequencies approaching the 
limit cycle frequency, is open to question. Second, there is a minimun 
amplitude of signal which can be detected corresponding either to the 
limit cycle amplitude or the zero-level width of the relay. Third, the 
bandwidth of the loop is limited by the sampling frequency, and signals 
above this frequency will not be passed; instead they appear as steps. 
(Ref. 19). 
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CHAPTER III 



DELTA-MODULATED INSTRUMENTS 



3. A General Comments 

The delta -modulated loops discussed in Chapter II are to be applied 
to two inertial instruments, the gyroscope and the accelerometer. 

The gyro to be used is the single-degree-of-freedom floated integrating 
gyro described by Draper, Wrigley and Grohe (Ref. 21). The 
accelerometer is a single -axis floated integrating pendulum type 
described by Hutchings and others (Refs. 4, 7, 8, 12, 17). 

3.A.1 Brief Description of the Instruments 

The integrating gyro, when behaving in an ideal manner, is 
described by the following differential equation: 



ja oa + 



CA„. = H W T . + M. 
OA sp IA tg 



(3-la) 



where J is the moment of inertia of the gyro float about its output 
axis; 

C is the damping coefficient of the fluid damper; 

H is the spin angular momentum of the gyro; 
sp 

Wja is the angular velocity of the gyro case about its input axis. 



with respect to inertial space; 

M. is the torque applied to the instrument float by the torque 
generator; 

and is the angle through which the float has turned with 

respect to the case avout the instrument output axis. 

A dot over a quantity (e.g. Aq^) denotes time differentiation. 




DELTA- MODULATED GYRO 

(a) 




DELTA-MODULATED ACCELEROMETER 

(b) 



Fig. 3-1 Delta-Modulated Instruments 
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Eq. (3 -la) may be rewritten as 




(3-lb) 



where 7 = J/C is the characteristic time or time constant of the gyro. 
The differential equation describing the accelerometer is 



where the symbols above are similarly defined; where P is the 
pendulosity of the accelerometer; and a^ is the non -gravitational 
acceleration of the instrument with respect to inertial space along 
its input axis. 

Eq. (3 -2a) may be rewritten as 



where again 7 = J/C. 

The torque generators are fed from a constant current source, with 
the current switched between two opposing windings to produce torque 
of different signs, or switched to a dummy load or off to produce zero 
torque. The required sign of the torque is determined at sampling 
times and is either opposite to the sign of the float angle or zero, 
depending on the logic of the loop. Therefore, the torque generator 
acts as a zero-order hold in the feedback path. 

3. A. 2 Form of the Loops 

The delta-modulated instruments may be represented as shown 
in Fig. 3-1. These loops are identical in form to each other and to 
the loops analyzed in Appendices A and B. The output pulses represent 
changes in the integral of the input, so that the gyro's output pulses 
indicate changes in the angular orientation of the gyro case about its 
input axis, and the accelerometer's output pulses indicated changes 
in velocity along the accelerometer's input axis. 

3.B Delta-Modulated Gyros 
3.B.1 Detailed Description 

A detailed physical description of the single -degree -of-freedom 



JA OA + CA OA Pa IA + M tg 



(3-2a) 




M 



(3 -2b) 
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floated integrating gyro is presented in Reference 21. Current 
practice is to center the gyro float in the case with a magnetic 
suspension system similar to that described in Reference 26, as 
well as to float it. Nominal values of the characteristics of some 
gyros are given in Reference 25. Fig. 3-2 shows a schematic 
drawing of the instrument. 

The differential equation governing the gyro may be derived 
from Euler's equations of motionfor arigid body about principal 
ax:es. It is assumed that the bearings of the gyro are infinitely stiff, 
so that only motions about the output axis need be considered. The 
bearing stiffness assumption is quite well satisfied for gyros even 
under high angular velocities. The set (IA, OA, SRA) forms a right- 
handed triad. The equation is 



M app = - H sp (W IA COS A OA " W 3RA sin A OA> + J < A OA + W ca> 



TA 3A n IA 

r T w 2 



TA SA 



M: 



‘ W SRA 2 > sin A OA cos A OA 




OA ' sln A OA ,W IA W SRA 


(3-3a) 


= M. - CA~. 
P tg OA 


(3-3b) 



where M is the torque applied to the instrument float about 
app 

its output axis; 

WgRA is the angular velocity of the gyro float about its spin 

reference axis, with respect to inertial space; 

W is the angular velocity of the case about its output axis, with 

Cel 

respect to inertial space; 

IjA is the moment of inertia of the gyro float about its input axis; 
IgA is the moment of inertia of the gyro float about its spin axis. 
These equations may be rewritten as 

7A OA + A OA " W 3RA sin A OA = “ P W IA COS A OA + 
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Fig. 3-2 Pictorial Diagram of Floated Single-Degree-of-Freedom 

Integrating Gyro 




M, 



- 7 W + 
ca 



C + (I SA 



x ia hw ia ' w sra ) sin a oa cos a oa 



/C 



+ (I SA • I IA ) (cos2 A OA ' sin2 A 



OA )W IA W SRA 



/C 



(3-3c) 

The last two terms on the right side of Eq. (3 -3c) are anisoinertia 
effects and are negligible in practical gyros. They will not be 
considered further here. 

Since A stays well below a degree in practice, the standard 
small angle approximations may be made in Eq. (3 -3c): 



7 ' A OA + A OA " — C^ W SRA A OA 
For W 



H 






IA - 7 W ca 



M, 

-c* 



(3-3d) 



‘'RA ~ ^ ca = 0* Eq. (3-3d) is identical with Eq. (3-lb). 



3.B.2 Rotational Coupling 

Eq. (3 -3d) shows that rotations of the gyros about all axes affect 
its motion. The coupling effects of these motions will be considered 
in detail. 

3.B.2.a Non -Reversing Rotations 

Examination of Eq. (3 -3d) shows that an angular velocity about 
the spin reference axis appears as an elastic restraint term in the 
gyro equation. This elastic restraint causes imperfect integration in 
the loop. As noted in Chapter II, Section 2.D.1.C, imperfect integration 
in the loop produces a dead zone, and inputs within this dead zone 
cannot be detected. This is a significant phenomenon in gimballess 
inertial reference equipment, since large angular velocities may 
occur about any axis. 

Angular acceleration of the case about its output axis produces 
a torque on the instrument float of the inertia, reaction category. As 
can be seen from Eq. (3 -3d), this torque is proportional to the gyro 
time constant and acts as a spurious input to the gyro. 

These torques may be compensated out either directly in the 
gyro or in the computer, if it is found that they seriously affect the 
performance of the gyro, in a manner similar to that used for 
compensation of compliance effects in present gyros. 
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3. B. 2. b Vibratory Rotations 

Vibratory rotations of the gyro about one axis at a time do not 
affect the performance of the instrument. When this type of motion 
occurs about two axes at a time, however, it can produce a drift. 

Consider sinusoidal motion of the case about the spin reference 
axis and the input axis. With sinusoidal motion, the gyro equation of 
motion, ignoring torque generator torque, is 



tA + A - Sf 

ta oa a oa ~cr 






(3-4) 

A solution for this equation using series techniques, which coverges 
for Wgj^ H g p/C < 1, has been given by Haff and Meltzer (Ref. 27, Ch. 5). 

A further discussion of the series solution and its convergence is given 
in Reference 29. 

Haff and Meltzer show that if n fgj^ = rnfj^, a constant drift is 



predicted by the series solution which they propose. A similar effect 
occurs if a vibration about the output axis occurs with the vibration about 
the spin reference axis. Further, if fgp^ is harmonically related to 



the clock frequency, drift may also be present, as can be seen by writing 
the Fourier series expansion for the torque generator torque. 

Angular vibrations of the gyro about any two axes produce a conical 
motion. Computation of the orientation of the gyro may be made in such 
a way as to introduce a constant angular velocity about the cone axis in 
the computed orientation, when none actually exists. Haff and Meltzer 
present a general treatment of this effect for a stabilized platform 
(Ref. 27, Ch. 3, Sec. 3. 4). It is clear from their discussion how the 
computation implied by a stable platform leads to this error. This 
coning error does not exist in gimballess systems. 

3.B.3 Linear Motions 

Many writers (Ref. 22, for example) have treated the effects of 
linear motions on gyros and they will not be considered here. 

3.B.4 Choice of Loop Parameters 

While the parameters of the loop must be chosen so that production 
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gyros and available electronic components may be used, a certain 
amount of latitude exists in choosing some parameters. 

The spin angular momentum of the gyro is chosen on the basis 
of the allowable drift rate. The higher the drift acceptable, the 
lower the angular momentum may be. 

The torque level, M, required in the torquer depends on the 
maximum angular velocity that the system must follow, W max . The 

product of angular momentum and W max gives the lower bound on M. 

For gyros with drift characteristics acceptable in inertial guidance 
applications, this bound may be quite high, and torquers other than 
microsyns may be required. 

The clock frequency may be determined from the maximum 
angular velocity expected and the increment size in the loop, A0, . 

The quotient W max /A0 gives the lower bound on clock frequency. 

If an ac signal generator is used, the carrier must be a multiple of 
the clock frequency. This requirement follows from the desirability of 
peak detection in the switch. If the carrier frequency is not a--multiple 
of the clock frequency, the phasing of the clock and the carrier will 
not remain constant, and the clock pulses cannot be adjusted to appear 
at the peaks of the carrier sinusoid. 

Since it is desirable to have the uncompensated mode of the loop 
as low as possible, or alternatively to have the relay zero -level width 
as narrow as possible, the gyro time constant, J/C, should be as 
low as possible. Since the polar moment of inertia of the float is 
usually fixed by the choice of angular momentum, the time constant 
depends only on the damping coefficient, C, which can be varied by 
varying the operating temperature of the instrument and the clearances 
between float and case. The higher the damping, the lower the time 
constant, but the smaller the angle through which the float travels. 
There is a certain value, A - n , below which the polarity of the signal 
cannot be determined, so that an indefinite increase in damping is not 
practicable. A short computation, suggested by J. E. Miller of the 
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M. I. T. Instrumentation Laboratory, gives an expression for the 
optimum damping in a loop with a two -level relay, considering only 
the dynamics of the gyro, which usually dominate. 

The differential equation governing the gyro is 



V CA OA = H sp W IA + M tg (3-5a) 

Integration of this expression over an integral number of sampling 
periods* and assumption of zero initial conditions for Aq^ gives 



ja oa + ca oa' h sp w ia 



- 8 . .,) 

ind 



(3-5b) 



where 0^ is the angle through which the gyro case has turned about 
its input axis; and 

0- nc j is the indicated value of 0^, obtained by counting output 
pulses. 

The quantity (0^ - ^ nc j) may be recognized as information stored 

in the float (see Ch. II, Sec. 2.C. 4). When the torque switches sign, 
the value of A is A m ^ n and the value of A is M(l-Or n )/(l+a' n )C. 

T 

s 

Here a = e 1 and the loop is assumed to be moding n-n. Thus at 
switching. 



j = JM(l-<* n ) + C 

S " ^p C1+ “ n> H sp Amin (3-5c) 

By differentiating Eq. (3 -5c) with respect to C, and setting the 

derivative equal to zero, an extreme value of I may be found, which 

s 

can be shown to be a minimum. The dependence of a and n on C may 
be neglected. 



or 




JM(l-a n ) + 

H dW) 
sp 



min 



H 



sp 



= 0 
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_ /(l-a n )JM 



(l+a n )A 



min 



(3-6) 



To show that this value of C indeed provides a minimum in I g , 



d 2 I 



dC 



— is computed: 



d 2 I 



q = 9 JM(l-a n ) >n 
! H 

sp 



dC' 



C 3 (l+a n ) 



which indicates a minimum. 

This optimum value of C has been chosen to minimize stored 
information. The corresponding time constant is 




r (l+<* n )JA 



min 



(l-cc n )M 



(3-7) 

For loops with a three -level relay, the damping should be chosen 
so that 

A = A 



min 

The expressions for A c are derived in Appendix B, Sec. B. 1. 



(3-8) 



ja oa + ca oa + kAqa - H sp W IA + M tg 



It is possible to show a relationship between the elastic restraint 
present in a gyro and the threshold of the instrument. The differential 
equation of the gyro is 

L \T A - TJ \\T 4 - 1 \/T 

(3-9) 

where K is the elastic restraint acting on the float from all causes, 
including electromagnetic reaction torque and the effect of angular 
velocity about the spin reference axis. If the inertia-reaction torque 
is neglected (it is usually much smaller than the damping torque) and 
the average angle of the float, Aq^, is considered, limit cycles may be 
ignored (See Ch. II, Sec. 2.D. 1. b) Under these assumptions, the solution 
of Eq. (3-9) is 

— h s P w ia -Hr 1 

A = S P z±L (\- r ^ ) 

OA K ' e ) (3-10) 
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In time T, the average angle of the float reaches a value A which 
represents one increment in the loop's output. This value is given by 



A = ^ r p 
A A 9 C s 

Writing Eq. (3-10) for t = T and solving for T yields 



T = 




A A9 

H W 

sp 



IA 



K) 



(3-11) 



T becomes very large as 



( A A0/W K approaches unity. 



For this 



term greater than or equal to one, T is not defined. Physically, this 
corresponds to the case where elastic restraint completely masks 
the input (See Ch. II, Sec. 2. D.l. c). Using only the first and second 
order terms in the series expansion of the logarithm, for small 
K» gives 



CA 



H sp W IA 



A 9 1 

(1 + T 



l A 9 



2 H W T A 
sp IA 



K) 



(3-12) 



Since MT 

s 




A 9, this may be written 



T 





A9 



W. 



IA 




(3-13) 



Thus the percentage error in time of occurrence of an output pulse 
is 



* =1 ^ j£ 

T 2 * WjA * C 



(3-14) 



By manipulating Eq. (3-14), many interesting relations between the 

time error and the loop parameters can be seen. The error in time 

decreases with elastic restraint and maximum input angular velocity. 

It is inversely proportiqnal to the actual input angular velocity and 

the square of the damping coefficient. Since the optimum damping 

coefficient depends inversely on A . , decreasing A . decreases 
r J min ° min 



6 



T* 
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3.C Delta -Modulated Accelerometers 
3.C.1 Detailed Description 

A detailed physical description of a pulse -restrained accelerometer 
is given in Reference 8. Present practice is to center the float in 
the case by a magnetic suspension system described in Reference 25. 
Reference 4 gives some parameter values for a typical accelerometer. 
Fig. 3-3 shows a schematic drawing of the instrument. 

The accelerometer's differential equation is 

j A oa + ca oa + Pa pRA sln A OA = Pa lA oos A OA " JW ca 

+ (I IA - I PA> W IA W PA + M tg (3-15a) 

where PA and PRA refer to the pendulum axis and pendulum reference 
axis respectively. (IA, PRA, OA) is a right-handed triad. 

The last term on the right hand side of Eq. (3 -15a) is an anisoinertia 
term which is negligible in practical accelerometers and it will not 
be considered further here. Making the small angle approximation 
for Aq^ and dividing by C gives 

P P . M, 



tA OA + A OA 



C a PRA A OA C a IA 



rW + 
ca 



tg 



(3 -15b) 

When a pj^ = W ca = 0, Eq. (3 -15b) becomes Eq. (3 -2b). As in the 
case of the gyro, infinite bearing stiffness has been assumed, and 
torques considered only about the output axis, 

3. C. 2 Rotational Coupling Effects 

As can be seen from Eq. (3 -15b), only rotations of the instrument 
about its output axis affect it, because of the assumed stiffness of 
the bearings. As in the case of the gyro, angular accelerations of 
the case about the output axis look like inputs. 

If the case motion about its output axis is vibratory and a vibratory 
linear motion along the pendulum reference axis occurs, an equation 
of the form of Eq. (3-4) governs the accelerometer motion. As noted 
above, if the frequencies of the two motions are harmonically related, 
a drift in the loop's output may appear. 
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Fxg. 3-3 Pictorial Diagram of Floated Single- Axis Pendulum 



3.C.3 Linear Coupling Effects 

Eq. (3-15b) shows that an acceleration along the pendulum reference 
axis appears as an elastic restraint. Imperfect integration results 
and a dead zone appears in the instrument. This is important in any 
application of this accelerometer since it may experience acceleration 
in any direction. 

If a vibration occurs along the pendulum reference axis in 
conjunction with a vibratory input, a drift may result as described 
above. A similar situation may arise if vibrations along the pendulum 
reference axis are related to the clock frequency. 

3.C.4 Choice of Loop Parameters 

Parameters for delta-modulated accelerometers may be chosen 
in a way similar to that used when choosing gyro parameters. Those 
which may be chosen are pendulosity, torque level, clock frequency, 
and damping. 

The pendulosity of the accelerometer is chosen on the basis of 
the level of torque uncertainties expected. The lower the torque 
uncertainty, the lower the pendulosity may be. This choice involves 
the same sort of compromise as the choice of angular momentum for 
a gyro. 

The required torque level is the product of the pendulosity and 
the highest level of input acceleration expected. This value is practically 
always lower than the corresponding value for a gyro. 

The clock rate is determined from the maximum acceleration 
expected and the velocity increments used. Their quotient gives the 
lower bound on clock frequency. 

The choice of damping coefficient for the accelerometer is made 
just as in the case of a gyro. For a limit cycling accelerometer, the 
optimum damping is 




(3-16) 



which gives a time constant 
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(3-17) 



For an accelerometer with a three-level relay, the damping coefficient 
should be chosen so that 



The remarks about elastic restraint made with respect to the 
gyro apply to the accelerometer as well. The equation corresponding 
to Eq. (3-14) is 



3.D Non-Ideal Components 

Up to this point, very little has been said about the non-ideal 
behavior of any components. The following sections consider some of 
the more serious hardware problems. 

3.D.1 Gyros 

In most applications, floated, single-degree -of -freedom integrating 
gyros do not experience significant angular velocities about any axis. 

In any environment in which they do, their performance will deteriorate, 
simply because they were not designed to take such punishment. A high 
angular velocity environment affects present-day gyros in two ways. 
First, they are not structurally strong enough to withstand the bending 
and torsion they experience. When torque is applied by the torque 
generator, the float experiences a torsional and bending stress because 
of the high spin angular momentum it has. In gimballess applications, 
the high torques involved may literally tear the instrument apart. 

A second effect of high angular velocities is high bearing loading. The 
bearings must be able to supply reaction forces strong enough to counter 
a torque equal to the torque generator torque. This cannot be done for 
long periods of time, but in a floated instrument the floatation fluid 
provides temporary additional bearing support, so that for short times 
the float will stay reasonably well centered. While at the present time 
floatation provides necessary additional support for the instrument , 




(3-18) 



€ T = 2 &ia C 



1 AV K 



(3-19) 
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magnetic supports are becoming capable of handling the entire load 
alone. 



3 .D.2 Torque Generator and Current Source 

In the analysis above, it has been assumed that the time constant 

of the torque generator was negligibly small. This assumption is 

violated as the torque from a given instrument is increased. Torque 

rises as the square of the input current. To increase this current, 

either the voltage must be increased or the impedance of the torquer 

decreased. Increasing the voltage is not an attractive solution, since, 

in general transistors drive the torquer, and their voltage limit is 

quite low. Lowering the impedance requires lowering the resistance, 

which in turn raises the time constant of the torquer, L/R. As the 

time constant of the torque motor increases, the effect of the lag 

makes itself felt in the moding plot. The value of the torquer time 

constant above which the lag must be considered is about T = '5'T . 

tg s 

More torque can, of course, be obtained by increasing the size of the 
torquer, but it is desirable to have small instruments. 



Microsyn torquers (Ref. 28) have been successfully used in pulse 
torquing applications. However, they are not high torque instruments, 
and their dependence on current squared makes them non-linear when 
used with a three -level relay. 

Flapper-type torquers (Ref. 23) provide about ten times as much 
torque as a corresponding size microsyn. Their main disadvantage 
is the presence of a fairly large negative elastic restraint, which 
makes the torquer positionally unstable. This might be eliminated by 
providing a positive mechanical or electrical spring, or using ac 
voltage for torquing and tuning out the restraint term. The ac torquer 
provides a double frequency term as well as a constant average, so 
that the frequency used must be well above the clock frequency, and 
above the mechanical break point of the system. 

Moving-coil, or d'Arsonval, instruments also provide large torques 
for small size. This type of torquer usually uses a permanent -magnet 
field. The difficulty with permanent magnets is that measurement 
techniques are not yet accurate enough to determine the stability of the 
field. 

A problem that all three kinds of torquers share is the stability 
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of their sensitivities. Good operation of the instruments requires that 
the same angular impulse be applied to the float for each electrical 
impulse applied to the torquer. If this is not the case, the instrument 
scale factor will vary with time, and errors will appear from this 
source. 

A related problem is the variation in torquer time constants with 
direction of torque. If the time constants are different, more impulse 
will be applied in the direction of the smaller time constant. Eventually 
a change in the output will be indicated although the input has been zero. 

To achieve the required accuracy in torquing . a very <■ '■ 

stable current source must be used. If a three -level relay is used, 
either the current must be turned off, or it must be fed into a dummy 
load. Turning off the current, while achieving great savings in power, 
causes stability problems in the current source above those already 
present. Ac current sources can be made more stable, but ac torquers 
are less sensitive than dc torquers. 

3.D.3 Signal Generator 

With ac signal generators, the output is a supressed-carrier wave. 
Thus, careful attention must be paid to the effects of circuitry which 
operates on the signal generator signal so that the gain for each 
sideband frequency is the same, and the phase shifts at each sideband 
are the negative of each other. If this is not the case, a quadrature 
component, which will cause a larger uncertainty in the position of the 
float with respect to the case, will be introduced into the signal. The 
end result will be an increase in the minimum angle of the float which 
can be detected. For this reason, the low-pass to band-pass trans- 
formations commonly used in ac compensation are not sufficiently 
accurate to be used in this application. 

As noted above (Sec. 3.B.4) when ac signal generators are used, 
the carrier and clock pulse phases are adjusted so that the clock pulses 
occur at the peaks of the carrier sinusoid. This adjustment is usually 
made by shifting the ac signal at either the primary or the secondary 
winding of the signal generator. For the reasons noted above this 
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phase shifting should be done on the primary to avoid adding a 
quadrature component to the signal generator and increasing the 
minimum discernable angle. 

Flapper-type (Ref. 24) and d'Arsonval signal generators provide 
attractive alternatives to microsyn signal generators, because of 
their small size. 
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CHAPTER IV 



KINEMATICS OF ROTATING COORDINATE FRAMES 



4. A Introduction 

In a conventional inertial guidance system, the inertial space reference' 
is preserved physically by the gimballed stable platform. The gimbals 
isolate the platform from the motion of the vehicle through servo drives, 
operating to null the output of gyros mounted on the stable platform. In 
a gimballess system, the inertial space reference may be preserved 
analytically in a transformation computer. This computer operates on 
a vector measured in body-fixed coordinates (the AV vector) to give the 
components of the vector in an inertially fixed frame of reference. The 
computer uses angular information obtained from gyros mounted on the 
vehicle to generate the rotational transformation. 

In this chapter, the properties of the transformation and the relations 
between the angular velocity of the body-fixed coordinate system with 
respect to inertial space and the rates of change of the transformation 
quantities will be adduced. 

There are essentially three methods of representing rotational trans- 
formations. Probably the most familiar is the matrix of direction cosines. 
This matrix, or alternatively, this second-order tensor, may have its 
elements written in terms of the nine direction cosines explicitly, in 
terms of a set of Euler angles, or in terms of the axis and magnitude of 
the rotation. Quaternions, invented by Sir William Hamilton in 1853 (Ref. 31), 
use four parameters to describe a rotation. The vector to be transformed 
is cast in quaternion form, pre -multiplied by a quaternion representing 
a rotation, and post-multiplied by the conjugate quaternion. The result 



is a quaternion representing the rotated vector. A variation on this, the 
Cayley -Klein parameters, uses the same four parameters in a slightly 
different form. The third transformation uses a vector to represent the 
rotation, and performs the transformation by means of purely vector and 
algebraic operations. 

These three types of transformation will be compared to determine 
which is most useful for machine computation. 

4. B Direction Cosi nes: Matrices and Tensors 
4.B.1 Direction Cosines 

The most familiar of the rotation transformations is the matrix 
transformation (Ref. 30; p. 95ff) 

vj = CV b (4-1) 

where Vj is a column vector with components measured in an 

inertial coordinate frame. 

V b is the same vector with components measured in a 
body-fixed coordinate frame, and 

C is the square matrix of direction cosines of the inertial 
axes relative to the body axes. 

The first component of Eq. 4-1 in vector notation is 



V I1 = % 



= (1 



II 



V,. 



H>i )v bi + (1 n 



W V b2 + (1 I1 



1 b3 )V b3 

(4-2) 



where ljj is the unit vector along the inertial X axis, etc. Since the 
vectors involved are unit vectors. 



% ' l bl = COS (1 I1 * W = C 11 (4-3) 

and in general, the elements of the matrix C are 

C ij = cos (1 I ' 1 bj ) = hi ' l bj (4-4) 

The rate of change of the elements of C may be found by differentiating 
Eq. (4-4) directly: 

cF C ij ill hi ' H>j) = ml 1 Ii, 



L . + 1 
bi 



Ii 
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or using a dot (• ) to denote time differentiation. 



C ij = C ij + 1 “j+2 " C ij + 2 W j+1 
.th 



(4-6) 



Note that w is the j component of angular velocity of the body -fixed 

J 

coordinate frame measured in body coordinates. 

Eq. (4-6) may also be derived by differentiating Eq. (4-1) directly: 

Ti = 6v b + cv b 



= C(V b + c' 1 ^) 



(4-7) 



-1 ’ — 

This is the Law of Coriolis in matrix form. If the term C CV b be 
identified with oo X V b , and the angular velocity matrix £2 be defined by 



0 



£2 = 



V 2 



GO, 



-CO, 



-GO- 



CO, 



then 



co X V, sA-£2V u 
b b 



c _1 cv. , 

b 



(4-8) 



(4-9) 



Thus, the differential equation for C is 

C = -C£2. (4-10) 

which expands to give Eq. (4-6). 

There are nine direction cosines. Since a rotation involves only 
three degrees of freedom, there must be six constraints on these nine 
quantities. These constraints arise from the fact that during a rotation, 
the magnitude of a vector remains unchanged. This condition requires 
that 



T -1 

C l = C 

T 

where C is the transpose of C 
and C 1 is the inverse of C. 

Another, more useful, form of the orthogonality condition is 



(4-11) 



C T C = CC T = I 



(4-12) 



where I is the identity matrix. Thus it is possible to write the nine 
direction cosines as functions of any three of them, using six independent 
orthogonality relations from Eq. (4-12). Another solution involves the 
use of Euler angles. 

4. B. 2 Euler Angles 

There are many ways of defining Euler angles; the definition 

shown in Figure 4-1 follows the most common usage (Ref. 30 p. 104ff. , 

and Ref. 39, p. 23, for example). The first rotation is about the Z axis 

of magnitude A^. The second rotation is about the displaced x axis, which 

becomes the line of nodes, by an amount A^. The final rotation is about 

the z axis by an amount A . In example being used, the capital letter 

z 

axes correspond to the inertially fixed frame, while the lower case axes 
are the body-fixed set. 



Z 



X 




LN 



Fig. 4-1 Euler Angles 

In terms of these Euler angles, the direction cosine matrix is 
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atn A l *in 



( 4 - 13 ) 



co« A_ co* A - sin A„ co* A, *in A 

Z Z 6 L Z 

-sin A^ co* A^ - sin Ay cos A^ cos A^ 
sin A z sin A^ 



- co* A z *in A z - sin A 7 cos A^ co* A z 

- sin A^ sin A^ + cos A^ cos A L co* Ay 

- cos Ay sin A^ 



sin A r cos A z 

'l 



'L' 
cos A t 



The rates of change of the Euler angles in terms of the 
angular velocity of the body-fixed set measured in the body-fixed 
coordinates may be determined by resolution, and are (Ref. 39; 
p. 24) 

. sin A cos A 

A z , z 

A r ? - co, — : — + CO ; 

Z 1 sm A^ 2 sin 

A, = w cos A + sin A. 

J — j L Z u Z 

A = oj cot A. sin A - w cot A T cos A + co 
zl L z2 L z3 

(4-14) 

This computation suffers from two disadvantages; trigono- 
metric functions must be determined, and a singularity exists 
at sin A^ = 0. The singularity corresponds to gimbal lock in 
a three-gimbal system (see Ref. 39, Figure 8-37, p. 25) and 
can be eliminated by including a fourth parameter in the 
description of the rotation, just as gimbal lock may be eliminated 
by the addition of a fourth gimbal. 

4.B.3 Tensors 

A third expression for the direction cosines is most 
conveniently cast in tensor form. It should be remarked here 
that since the transformation involves only Cartesian coordinates, 
no distinction between covariant, contravariant, and physical 
components of vectors and tensors need by made (Ref. 35, p. 170). 
The direction cosine tensor is 



C. . = cos A 6. . + (1- 
ij iJ 



cos A) r.r. - sin A e... r, 
i J ijk k 



(4-15) 
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where the rotation which carries the body -fixed axes into the inertial 
axes is about an axis with direction cosines r^ (the same with respect 
to either coordinate system) with magnitude A. Here is the 
Kronecker delta 



6..=/ lifi = j 

1J \0 if i f j 



(4-16) 



e ijk * s c y clic tensor 



e ijk 



+1 if ijk positive cyclic (e. g. 123) 

= {-l if ijk negative cyclic 

0 if i = j, j = k, or i = k (4-17) 

and the convention of summing on repeated indices is observed (Ref. 34, 
p. 96 for example). Using the tensor formulation, Eq. (4-1) is 



V Ii C ij V bj (4-18) 

and Eq. (4-10) becomes 

^ij ^ik^jm^m (4-19) 

4. C Four Parameter Representation 

The direction cosine computations are six -fold redundant; the Euler 
angle computations carry no redundancy but have a singularity. The use 
of four parameters to describe a rotation eliminates the singularity and 
diminishes the redundancy. There are two practically identical ways of 
describing a rotation with four parameters. These are Hamilton's 
quaternions (Refs. 31-33, 36, 38, 40) and Cayley-Klein parameters 
(Refs. 30, 40). 

4.C.1 Quaternions 

Hamilton developed quaternions as an extension of complex numbers. 
In representing a vector in two dimensions as a complex number, the 
unit real number 1 represents one direction, and the unit complex number 
i = n/- 1, the other. In analogy with this, Hamilton let the unit space 
vectors be i, j, k, with multiplication rules similiar to those for si- 1: 

. i ^ = y? = -1 ij = -ji = k, etc. (4-20) 

Thus, a space vector may be represented as V = iV^ + jV2 + kV^, and 
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a general quaternion, with the unit real 1 the fourth direction, as 

q = a + ib + jc + kd (4-21) 

Conjugation is defined by 

q' = a - ib - jc - kd (4-22) 

and length is given by 

i 2 = q 5 q = a 2 + b 2 + c 2 + d 2 (4-23) 

A simpler notation follows from representing a quaternion as a scalar 
and a vector: 



q = (X, p) = X+ ip x + jp 2 + kp 3 

Thus; > 

q' = (X, -p) 

and if q^ = (X^ Pj) and q 2 = (X 2 , p 2 ), then 

q = (q^2^ = ^ 2 *^ 2 ) 

(XjX 2 - Pj.P 2 ‘ *1^2 + ^ 2^1 + ^ 1 X ^ 2 ^ 



(4-24) 

(4-25) 



(4-26) 



Note that in general multiplication is not commutative, . Two special 
quaternions may be defined, the unit quaternion q^, and the zero quaternion 

V 



q u = U. o) q o = (0,6) 



(4-27) 



Obviously 



qq u = q u q - q q G q = qq Q = q o u-28> 

Further, if q f q Q , there exists a unique inverse or reciprocal q ^ such 

+ -1 "I 

that qq = q q = q u = 

-1 1 * 

q = — o q 

I ql (4-29) 

Two relations follow from the definitions of conjugate and inverse: 

(q x q 2 ) * = q 2 *qi* ^1^2* = q l q 2 (4-30) 

If the quaternions qj and q fe ar{ , de£ined by 
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“i = <a r V % ° (a b- "V 



x2 , 2 _ , 

A + p - 1 



and related through the quaternion q having unit magnitude. 

q = (X,^) 

by the transformation 

qi = *%* 

then 



(4-31) 

(4-32) 

(4-33) 



(a r Vj) = [ a b , V b + 2 Xp X V b + 2p X.(p X V^)] (4 _ 34) 

The scalars of qj and q b are equal, and the vectors are related by 

V T b + 2Xp X V b + 2p X (p X— b )., (4 35) 

Writing out Eq. (4-18) in vector notation, and using r as the unit vector 
along the axis of rotation gives 

Vj = V b cos A + (r- V b ) r (1-cos A) + r X V b sin A 

* V b + r X V b sin A + [ (r* V b )r - (r* r)V b ] (1-cos A) 

= V b + r X V h sin A + r X (r X V b )(l-cos A) 



V b + 2 sin ^ cos ^ r X V b + 2 sin 2 r X (r X V b ) 



A . 



(4-36) 



If X in Eq. (4-35) is identified with cos -w- in Eq. (4-36), and p with 
- A z 

r sin-^-, then the transformation Eq. (4-32) corresponds to the rotation 
of T b about an axis r through an angle A. Writing q out gives 

A - A A A 

q = (cos ~ 2 , r sin-^) = cos-?, + sin-^ ( ir i + J r 2 + kr V 

Extending Euler's formula for exponentials of complex quantities 
suggests the notation 

A 

~o ( ir 1 + J r 2 + kr 3 ) A A 

e ^ = cos-4- + sin -4 (ir. + jr„ + krj 

£ £ l £ 6 (4-38) 

2 2 2 

with r^ +r 3 + r 3 = 1> so that the transformation Eq. (4-33) may be 

written 
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i (ir l + J r 2 + kr 3 ) " T (ir l + jr 2 + kr 3 ) 

q x = e q fe e 

The quaternions representing successive rotations may be combined as 
follows: 

211 2 21 21 (4-40) 

The non- commutativity of finite rotations is evident from the non- 
commutativity of quaternion multiplication. 

The rates of change of the rotation quaternion parameters must now 
be related to the angular velocity of the body-fixed coordinate system 
with respect to inertial space. Let q be the quaternion representing the 
rotation that carries the body-fixed space into the inertial frame at time 
t. The quaternion representing the opposite rotation is then q . Let 
the quaternion representing the rotation carrying the body-fixed frame into 
inertial space at time|t + A t|be q 1 . The quaternion representing the 
rotation of the body-fixed axes during the time At is therefore q'q*: 

AA AA - 

q'q* = (cos — > sin j r^.) = (X + AX, p+ A p)(X,-p) 

= [ (X, p) + (AX, Ap)] (X, -p) 

= (1,0) + (AX, A p)(X,-p) 

(4-41) 

or 

(cos 1, sin ^J“ r At^ = (AX, Ap)(x, -p) (4-42) 

Dividing this equation by At and allowing At to approach zero gives 



lim 1 
At -*0 2 



cos 



AA 



1 



(- 



A A/2 



AA, - 
~Tl)r 



sin 



AA 



AA 



1 lim 

2 (0,, r 



At, At AA/2 At 
AA 






t-0 At At 



-) = (X,p )(X, -p) 



(4-43) 



Since ^ At ^ Ec l- (4-43) becomes 
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2 (0,u>) = (X, p MX, -p) 



(4-44) 



or, post -multiplying by (X,p) 

(X, p ) =2 (0,w )(X,p) 

whence 




and 




(Xw + wXp) 



(4-45) 



(4-46) 



4. C.2 Cayley-Klein Parameters 

The foregoing discussion provides a method for describing a 
rotation with four parameters. One is redundant; the redundancy is 
removed by the constraint that the quaternion q = (X,p) have unit 
magnitude [ Eq. (4-32)] . This formulation suffers from the fact that 
quaternions are not at all widely known. To obviate this difficulty, 
Cayley-Klein parameters may be introduced; they retain the simplicity 
of quaternions while eliminating Hamilton's i, j, k. Defining the Cayley- 
Klein parameters by 

a = X " iP 3 

& = -P 2 - 1 P 3 (4-47) 

T - +P 2 - ip 3 
b = X + ip 3 

and from them, forming the complex matrix Q, 





Q = 



gives a matrix which is unitary and has determinent + 1 : 
*T 

Q = Q - Q 1 



I Ql 



<*** + / 5 / 3 * = X 2 + Po + p 2 + p? = X 2 + p 2 = +1 



(4-48) 



(4-49) 



’3 K 2 K 1 

❖ T 

Here, Q is the complex conjugate of Q, Q is the transpose of Q,and 

A 

Q is the adjoint of Q. When using Cayley-Klein parameters, vectors 
are represented as Hermition, or self-adjoint, matrices: 
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-« » t i A f = 



V 3 

IVl + iv 2 




(4-50) 



Note that the trace of the matrix V is zero, and the determinant of V is 
the negative of the square of the magnitude of the vector V: 

I VI = -(V? + V 2 + V 2 ) - -V 2 (4-51) 

Using the matrix Q, a similarity transformation may be performed 
on the matrix yielding another matrix 

Vj = QV b Q A (4-52) 

Under similarity transformations, the Hermition property of a matrix is 
preserved, as are the trace and determinant. Therefore, Vj must be of 
the same form as [ Eq. (4-50)] , and further 



- V b= IV b l - IVjl "Vj 



(4-53) 



Eq. (4-53) is the orthogonality condition. Thus the transformation 
Eq. (4-52) is equivalent to the transformation Eq. (4-33). 

The differential equations for the Cayley-Klein parameters follow 
directly from the defining equation [ Eq. (4-47)] and the rates of change 
of the quaternion parameters [Eq. (4-46),] , 

4.D Vector Representation 

Another useful representation of a rotation is the vector cr 

cr = r tan A (4-54) 

where r and A are defined above. In terms of the scalar and vector of 
the rotation quaternion 



a = r 



. A 
sm -4 

A~ 

cos -4 



= r 



• A , 
sin T 

O 2 A 
2cos 



1 x A 

1 + cos ~2 



1 + X 



(4-55) 



The inverse relations are 



. A 0 2 A . . 

X = cos — = 2 cos — -1 = 



sec 



2 A 



-1 = 



2-sec 2 4- 

4 

2 A 
sec — r 
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and 



(4-56) 



2-<l + tan 2 ^) 1 _ (J 2) 

, . , 2~K . , 2 

1 + tan 1 + o 



p = (1 + X) a = 



1 + a 



Substituting these expressions for X and p into Eq. (4-35) gives 

2 

“ ( — ~— ■ ) XV, + 2 ( 

1 + a 
or 



v i= v b + 



2 ( T^T )xv b + 2 <774* X I ( T 3 - 2 ) x v bl 



1 + o 



vl: /.; 

r\ 



8 



T T =T h + 4(1 . : , .^ axv h +■ 

1 b (1 + a ) b (1 + a *) 



n aX(ffXV b ) (4-57) 



The vector representation of a rotation requires only three numbers; 
there is no redundancy. There is a singularity, however, at|A] = 2ir , 
corresponding to a full rotation. This singularity may be eliminated by 
defining another vector o': 




(4-58) 

o' = i V - = ( - r) cot ^-= (-r) tan (y - 'y) 

tan 



o' = (- r ) tan ■ 2n ~ (4-59) 

Since a rotation of A in one direction is equivalent to a rotation of 
2n -A in the opposite direction, a and cr' describe the same rotation, and 
Eq. (4-57) holds with o' substituted for a. 

The differential equation relating the rate of change of o to the 
angular velocity of the body -fixed coordinate frame may be obtained by 
differentiation of Eq. (4-55): 

5 = _i ip 

1 + x a + x :) 2 



70 



_ 1^ . Aco + coX p ( w • p) p . 

'2 l 1 + X (l + rf 1 

1 1 -o 2 

= sr [ — = — co + to X a + (w ct)ct] 

* (4-60) 

1 l 2 

= J [ w + CO X a + (coxa) X a + aw] 

y 2 _ 

o' = 2"[ * 0 | co + to X a + (co Xo) X a] 

Clearly, the same differential equation holds for o': 

• 2 

o' = j [■— ' J> + wX o' + (coXa 1 ) X o'] (4-61) 

4. E Comparison 

The preceding three sections of this chapter have presented several 
methods for transforming a vector from one Cartesian coordinate frame 
to another, and derived the differential equations for the transformation 
parameters in terms of the angular velocity of one frame with respect 
to another. These transformations and their differential equations are 
collected below for comparison. 

4. B Direction Cosines 
Transformation 

Vj = CV b (4-1) 

Direction Cosine Differential Equation 

C = -Cfl (4-10) 

Six -fold Redundant, with Six Constraints 

C T C = I (4-12) 

Euler Angle Differential Equation 

. sin A cos A 

A Z ~ W 1 sin A^ + ^2 sinA^ 

A t = co cos A + to sin A (4-14) 

LI z 2 z 

A = -co cot A t sin A cot A T cos A + co 

z 1 L z2 L z3 
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No redundancy; singularity at sin = 0 

4»C Four Parameters 

Quaternion Transformation 

V~ = Vfo + 2 Xp X V h + 2 p X (p X v b ) 

A A 

X = cos ~2 p = r sin 



Cayley -Klein Transformation 



”11 "12 
- V I 3 




-p, - ip, 
X *ip 3 




b2 



X + ip 3 P 2 * lPj] 
?2 + Vi X - ip, | 



(4-35) 

(4-37) 



( 4 - 52 ) 



Differential Equation 

' 1 - - 
X - - P 



1 



p - y (Xu) +ooXp ) 

Once redundant with one constraint 



X 2 + p 2 = 1 



4.D Vectors 



Transformation 



(4-46) 



V, = V, + 4 1 " ° 0 o X V, + — - 

1 b (1 + o 2 ) 2 b (1 + o “) 



2 s 2 



(4-32)(4-49) 



o X(ct XV b ) 



(4-57) 



o = r tan 



(4-54) 



or 



V I = V b + 4 



l-o 



,2 



(1 + o' ) 



— 5 - 9 -a X V. + 
,| 2\2 b 



— s-j-g-S' x (5- x v b) 

(l + o'T ' 



o = 



(4-57a) 

(4-58) 
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Differential Equation 




2 



co + co X a + (ojXct) X ct] 



( 4 - 60 ) 



.2 



a 1 = — co + coXj' + (coXa')X a' ] 



( 4 - 61 ) 



These transformations will now be examined with a view to determining 
which is most suitable for machine computation. The basic criterion to 
be used in this determination is the size of computer required. Further, 
since the information being dealt with is incremental, incremental 
computation will be the basis of decision unless general propose computation 
seems particularly appropriate. 

The most straightforward computation employs the direct direction 
cosine formulation. The transformation requires nine incremental 
integrators, since each component of the vector to be transformed has 
magnitude AV. The generation of the direction cosines requires eighteen 
incremental integrators. If Euler angles are generated, then the sine 
and cosine of each angle is also required; further, the sines and cosines 
must be combined to give the direction cosines as well as the factors in 
the differential equations. The use of incremental devices requires 
forty-six integrators. This does not handle any problems raised by 
the singularity at sin = 0 

Twelve incremental integrators are necessary to generate the four 
parameters used in the quaternion and Cayley-Klein parameter representation. 
If the direction cosines are formed from these parameters, an additional 
twelve integrators are required. Nine more complete the transformation. 
Thus, thirty-three incremental integrators are required in all, as opposed 
to twenty -seven for the direction cosine computation. If the formal 
Cayley-Klein similarity transformation is used, the first multiplication 
can be performed by twelve incremental integrators. An arithemetic 
computation seems best for the second matrix multiplication; it requires 
twelve multiplications and nine additions. 

Twenty-one incremental integrators are required to generate the 
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three components of the vector used to represent a rotation; the vector 
transformation requires twenty-five more. In addition, an arithmetic 
section is required to perform the branching operation used to remove 
the singularity. 

The foregoing discussion indicates that the direction cosine 
computation using p. total of twenty-seven integrators is most compact. 
For this reason, this formulation will be the only one considered in 
designing the transformation computer. 
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CHAPTER V 



TRANSFORMATION COMPUTER 



5. A Introduction 

In Chapter IV, the various ways of representing coordinate ro- 
tations were derived and compared to determine their suitability for 
machine computation. As a result of this comparison, the matrix trans 
formation was chosen to be instrumented. The rotation matrix is to be 
generated by direct computation of the direction cosines. The matrix 
differential equation to be instrumented is 



where C is the matrix of direction cosines, a dot (• ) over a quantity 
denotes time differentiation, and is the angular velocity matrix. 



The angular velocity is resolved into body axes. 

Eq. (5-1) represents nine separate equations; there are three 1 
sets of the form 



C = -CO 



(5-1) 



— — uj ^ 





(5-2) 



c il - C i 2 w 3 - C i 3 w 2 
C i2 = C i3 w l " C il w 3 
C i3 = C il w 2 ' C i2 w l 



(5-3) 



Multiplication of Eq. -(5-3) by dt and integration of the result yields 



5.B Algorithms Available 

Equations (5-4) show that each direction cosine is the difference 
of two integrals, each the integral of a direction cosine with respect 
to rotation about a body axis. Since the gyro outputs are equal incre- 
ments of rotation about body axes, the integrals may be approximated 
by quadrature rules using equal A0 increments. 

5.B.1 Quadrature Rules 

Almost any standard text on numerical analysis (Refs. 41, 

46-48, 52, 54, for example) gives the more common quadrature 
rules. Since only the integrand will be available, no rules employing 
derivatives of the integrand, such as Taylor' s series, will be con- 
sidered. The most common quadrature rules using present and past 
values of the integrand are the rectangular rule. 



C il ” f C i2 d0 3 * J C i3 d0 2 
C i2 = J C i3 d0 l - / C il d0 3 
C i3 = / C il d0 2 ~ f C i2 d0 l 



(5-4) 



0(n) 



Jc(0)d0 = A0C[0(n-l)] 
0(n-l) 



(5-5) 



the trapezoidal rule. 



9(n) 

j , C(0)d0=^-{;^9(n-l)]+C[0(n)] } 
0(n-l) 



(5-6) 



and Simpson' s one-third rule. 
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